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Abstract 

Li1+jCFeP04  (0  <  x  <  3)  as  anode  material  for  lithium  ion  batteries  has  been  studied  using  ab  initio  calculations.  Results  show  that  large  amount 
of  lithium  ions  can  be  intercalated  into  LiFeP04  host.  The  structure  changes  continuously  when  the  first  two  Moles  of  lithium  ions  (. x  <  2)  are 
intercalated  into  the  LiFeP04  host,  accompanied  by  large  volume  expansion  (37.4%  and  25.4%  for  the  first  and  second  Mole).  The  final  product 
of  Li3FeP04  possesses  a  stable  chained  structure,  which  is  favorable  for  storing  even  more  lithium.  In  the  same  time,  lithium  ion  diffuses  in  a 
three-dimension  pathway  within  the  chained  structure.  The  unit  cell  volume  increases  only  by  4.9%  from  Li3FeP04  to  Li4FeP04,  and  the  chained 
structure  keeps  unchanged. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  phospho-olivine  LiFePCL  has  been  extensively  studied 
as  cathode  material  for  lithium  ion  batteries  since  the  pioneering 
work  by  Goodenough  and  co-workers  [1].  With  advantages  like 
high-energy  density,  low  cost  and  good  environmental  compati¬ 
bility,  it  has  attracted  great  attention  from  the  standpoint  of  both 
academy  and  industry. 

However,  the  intrinsic  low  electronic  and  ionic  conductivity 
of  this  material,  which  results  in  rather  poor  rate  performance 
during  cycling,  is  the  major  problem  that  inhibits  its  wide  appli¬ 
cation  in  lithium  ion  batteries  industry.  Many  studies  have  been 
done  to  improve  the  dynamic  properties.  It  is  shown  that  the 
electronic  conductivity  of  LiFePCL  could  be  largely  improved 
through  Li  site  doping  [2,3],  but  the  dopant  cations  at  the  Li 
sites  might  block  the  lithium  one-dimensional  diffusion  pathway 
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[4,5].  Numerous  basic  researches  have  also  been  carried  out  to 
understand  both  electron  transport  mechanism  and  lithium  ion 
transport  process  during  cycling  [6,7].  Simultaneously,  various 
studies  have  also  been  concentrated  to  improve  the  rate  perfor¬ 
mance  via  decreasing  the  particle  size  or  improving  electronic 
contact  [8,9].  Although  many  efforts  have  devoted  to  improve 
the  rate  performance,  it  is  still  an  open  difficulty  to  both  academic 
and  industrial. 

Recently,  Kalaiselvi  et  al.  [10]  have  successfully  shown 
that  LiFePCL  compound  can  be  employed  as  anode  material 
for  rechargeable  lithium  batteries.  The  reversible  capacity  of 
about  300  mAh  g-1  (after  20  cycles  of  charge/discharge)  of 
LiFePCL  anode  indicates  that  ~2  lithium  ions  can  be  reversibly 
inserted  into  one  LiFePCL  formula,  and  becomes  LiFePCL. 
Later,  the  capacity  of  LiFePCL  as  anode  material  is  further 
improved  by  doping  with  small  amount  of  Cu  or  Sn  at  Fe  sites 
[11].  The  cycling  capacity  of  the  doped  LiFePCL  as  anode 
material  is  enhanced  to  ~400mAhg_1,  indicating  that  even 
more  lithium  ions  can  be  reversibly  stored  in  Lii+;cFePCL  host. 
Although  it  has  been  shown  that  LiFePCL  as  anode  possesses 
very  good  capacity  and  cycling  performance,  the  lithium 
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storage  mechanism  behind  the  performance  is  not  identified 
by  the  authors.  For  example,  structural  transition  and  change 
of  Fe  valence  state  during  lithium  intercalation  process  are  not 
studied,  which  are  obviously  important  for  understanding  of 
the  lithium  storage  mechanism. 

In  the  present  work,  we  are  going  to  investigate  the  basic 
physical  and  chemical  problems  of  LiFePCU  as  anode  material 
for  lithium  ion  batteries  from  theoretical  aspect  of  view,  namely 
from  density  functional  theory  (DFT)-based  ab  initio  calcula¬ 
tions.  The  structural  and  electronic  properties,  charge  transfer, 
lithium  diffusion  upon  lithium  intercalation  are  studied  in  detail. 

2.  Computational  details 

We  have  performed  the  DFT-based  ab  initio  calculations 
using  DACAPO  code  [12].  This  code  solves  the  Kohn-Sham 
equations  within  the  pseudopotential  approximation  whereby 
the  electrons  are  described  in  the  generalized  gradient  approx¬ 
imation  (GGA)  [13].  The  valence  and  core  interactions  are 
described  by  the  Vanderbilt  pseudopotentials  [14].  The  valence 
electrons  wave  functions  and  the  augmented  electron  density 
are  expanded  in  a  plane  wave  basis  sets  with  cutoff  energies 
of  25  and  140  Ry,  respectively.  The  Perdew-Wang  exchange- 
correlation  functional  (PW91)  [15]  is  used  for  the  calculation  of 
the  electron  exchange  correlation  energy. 

LiFePC>4  has  ordered  olivine  structure  (space  group  Pnma ), 
with  lattice  constants  of  a  -  10.334  A,b  =  6.0105  A,  c  =  4.732  A 
[1].  The  unit  cell  contains  4  LiFePCU  formulas,  including  28 
atoms  in  all.  For  structural  relaxation  and  total  energy  calcula¬ 
tions,  Monkhorst-Pack  scheme  [16]  3  x5  x7  k points  sampling 
is  used  for  the  integration  in  the  irreducible  Brillouin  zone. 

The  migration  of  lithium  ions  are  simulated  with  nudged  elas¬ 
tic  band  (NEB)  method  [17,18].  The  migration  energy  barriers 
are  calculated  within  a  dilute  vacancy  situation,  with  only  one  Li 
vacancy  and  otherwise  all  occupied  Li  atoms  in  the  unit  cell.  As 


Table  1 

Lattice  constant  (in  A),  unit  cell  volume  V  (in  A3),  magnetic  moment  of  Fe  atom 
MOM  (in  /i b),  and  total  energy  difference  A E  (in  eV  per  unit  cell)  of  FeP04 
and  LiFeP04 


a 

b 

c 

V 

MOM 

A  E 

FeP04 

Non- spin 

9.719 

5.731 

4.720 

262.90 

- 

Reference 

FM 

10.009 

5.980 

4.856 

290.65 

4.02 

-4.77 

AFM 

10.022 

6.021 

4.850 

292.67 

3.93 

-5.53 

LiFeP04 

Non-spin 

10.290 

5.982 

4.673 

287.65 

- 

Reference 

FM 

10.651 

6.206 

4.817 

318.36 

3.55 

-5.22 

AFM 

10.627 

6.187 

4.800 

315.61 

3.51 

-5.38 

A E  is  calculated  with  non- spin-polarized  total  energy  as  reference. 


the  NEB  calculations  are  computational  expensive,  only  gamma 
point  is  used  for  the  NEB  calculations. 

3.  Crystal  structure  upon  lithium  intercalation 

In  the  present  work,  we  first  optimized  the  theoretical  equi¬ 
librium  lattice  constant  of  FePCL  and  LiFeP04.  The  results  are 
listed  in  Table  1 .  The  theoretical  lattice  constants  for  the  spin- 
polarized  calculations  are  larger  than  the  non- spin-polarized 
one.  From  comparing  the  total  energy  of  the  system  under  dif¬ 
ferent  magnetic  configurations,  we  conclude  that  the  bulk  spin 
configuration  is  antiferromagnetic  (AFM),  in  agreement  with 
the  results  in  Ref.  [19].  For  the  AFM  configuration,  the  unit  cell 
volume  of  FeP04  increases  by  about  7.8%  after  lithium  ions  are 
intercalated  to  become  LiFePCL. 

Fig.  la  shows  schematically  the  crystal  structure  of  LiFePCL. 
Within  the  unit  cell,  four  Li  atoms  take  the  4 a  positions  while  Fe 
and  P  atoms  occupy  part  of  the  Ac  positions  and  O  atoms  occupy 
part  of  Ac  and  8 d  positions  [1,20].  From  the  crystal  structure,  we 
can  see  that  there  are  a  lot  of  vacant  sites  available  for  lithium 
storage,  including  Ab,  Ac,  and  8 d  sites.  When  more  lithium  atoms 


Li  4aFeP04  Li  4tfLi4/,FeP04  Li  4aLi4TeP04 

Fig.  1.  Schematic  view  of  the  atomic  structures  of  (a)  LiFeP04,  (b)  Li4aLi4/7FeP04,  and  (c)  Li4aLi4CFeP04.  The  largest  red  spheres  are  the  O  atoms,  the  middle 
sized  spheres  indicate  the  Fe  (grey)  and  P  (green)  atoms,  and  the  smallest  spheres  are  lithium  atoms,  with  purples  color  at  the  4 a  sites,  white  at  4b  sites,  and  grey  at 
4c  sites.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Table  2 

Lattice  constant  (in  A),  unit  cell  volume  V (in  A3),  magnetic  moment  of  Fe  atom 
MOM  (in  /xb),  and  total  energy  difference  A E  (in  eV  per  unit  cell)  of  Li2FeP04 
optimized  for  Li  occupy  4c  and  4b  sites 


a 

b 

c 

V 

MOM 

A  E 

Li4aLi4/?FeP04 

Non-spin 

12.027 

6.430 

4.923 

380.69 

- 

Reference 

FM 

12.194 

6.755 

5.077 

418.20 

3.79 

-8.94 

AFM 

12.340 

6.966 

5.118 

439.94 

3.75 

-9.78 

Li4aLi4CFePC>4 

Non-spin 

11.905 

6.558 

5.135 

400.93 

- 

-0.17 

FM 

11.960 

6.890 

5.246 

432.31 

2.99 

-11.41 

AFM 

11.932 

6.920 

5.252 

433.66 

3.03 

-11.71 

A E  is  calculated  with  non- spin-polarized  Li4flLi4^FeP04  as  reference. 


are  inserted  into  the  LiFePC>4  lattice,  we  first  need  to  know  which 
site  will  be  taken  by  these  lithium  atoms.  4 b  or  4c  sites  are  the  first 
choice.  From  the  symmetry  of  crystal,  the  fractional  coordinate 
of  one  4b  site  is  (0,  0,  0.5).  However,  for  the  4c  sites,  only  the 
y-coordinate  is  fixed  to  be  0.25,  the  x-  and  z-coordinates  are 
flexible. 

To  simulate  the  structure  of  L^FePCU,  which  is  formed  by 
intercalation  additional  four  lithium  atoms  in  one  LiFePC>4  unit 
cell,  we  try  to  place  these  additional  lithium  atoms  into  either  4b 
or  4c  sites,  and  the  unit  cells  are  denoted  as  Lk^LL^FePC^  and 
Li4aLi4cFePC>4,  respectively.  The  calculated  results  are  listed  in 
Table  2.  From  the  calculated  total  energy,  we  find  out  that  the 
optimized  bulk  magnetic  configuration  of  both  LL^Lf^FePCU 
and  Li4aLi4cFeP04  are  antiferromagnetic.  The  total  energy  of 
the  unit  cell  for  lithium  atoms  at  the  4c  sites  is  about  2eV 
lower  than  that  with  lithium  taking  the  4b  sites.  This  indi¬ 
cates  that  the  4c  sites  are  more  favorable  for  lithium  storage 
at  the  initial  stage  of  lithium  intercalation  process.  However, 
in  both  cases  we  observed  a  very  large  lattice  expansion  from 
LiFePCU  to  Li2FeP04.  For  the  energetically  most  favorable 
AFM  Li4aLi4cFeP04  configuration,  the  volume  of  the  unit  cell 
increases  about  37.4%,  from  315.61  A3  to  433.66  A3.  This  is  not 
favorable  to  the  stability  of  the  material  during  the  intercalation 
process. 


Presented  in  Fig.  lb  and  c  are  the  optimized  atomic  struc¬ 
tures  of  LiFePCU  with  AFM  spin  configuration.  In  the  case  of 
LiFeP04,  the  FeC>6  octahedra  and  PO4  tetrahedra  are  connected 
with  shared  edges  and  comers,  and  the  crystal  is  compactly 
organized  in  three-dimensions  [1],  which  ensure  the  structure 
to  be  very  stable.  However,  in  the  case  of  L^FePCU,  either 
LL^Lk^FePCU  or  LL^LLtcFePCU,  the  structure  is  quite  loose, 
and  the  unit  cell  volume  is  much  larger  than  that  of  LiFePC>4, 
as  discussed  above.  In  LL^Lf^FePC^,  FeC>6  octahedra  struc¬ 
tures  are  transformed  into  tetrahedra,  coordinated  with  one  P 
and  three  O  atoms.  PO4  tetrahedra  structures  are  also  changed, 
with  P  and  Fe  atoms  directly  coordinated.  Due  to  these  changes, 
the  unit  cell  is  organized  in  two  rings.  Each  ring  is  formed  with 
two  P-tetrahedra  and  two  Fe-tetrahedra.  The  structure  is  formed 
with  these  cluster  rings.  For  each  ring  structure,  the  unshared 
(one-coordinated)  O  atoms  are  stretched  into  lithium  layers, 
providing  certain  electron  affinity  and  stabilizing  those  lithium 
atoms.  In  Li4aLi4cFeP04,  the  FeC>6  octahedral  structures  also 
disappeared.  The  PO4  tetrahedra  structures  are  chained  along 
^-direction.  Those  tetrahedra  are  connected  to  each  other  with 
Fe  atoms,  which  are  bonded  by  three  O  atoms  and  the  FeC>3 
structure  is  arranged  within  a  plane. 

In  both  LL^LL^FePCU  and  LL^LL^FePCU,  the  crystal  struc¬ 
ture  is  substantially  changed.  The  unit  cell  volume  is  largely 
increased  and  the  structure  becomes  relatively  loosen  and  unsta¬ 
ble.  As  for  electrode  material  for  lithium  batteries,  this  is 
obviously  unfavorable  for  the  cycling  performance.  However, 
the  loose  structure  provides  large  space  for  lithium  storage, 
which  enables  the  high  capacity  as  electrode  material. 

Fig.  2  shows  the  optimized  atomic  stmctures  of  Li3FeP04 
with  AFM  spin  ordering,  which  is  also  the  energetically  most 
favorable  configuration,  as  shown  in  Table  3.  By  contrast 
with  the  structures  of  LiFePCU  and  L^FePC^,  the  structure 
of  the  material  and  symmetry  of  the  atoms  are  completely 
different.  The  structure  of  L^FePCU  is  constructed  with  quasi- 
one-dimensional  chains  along  the  z-axis  direction.  The  chain 
structure  is  formed  in  a  zigzag  way  by  alternative  Fe  and  O 
atoms,  with  one  PO3  group  bonded  to  each  Fe  atom  along 
the  chain.  Within  the  chained  structure,  the  atoms  are  orga- 


Fig.  2.  Atomic  structures  of  LbFePCL  (left)  and  the  chained  structures  in  the  crystal  viewing  from  different  direction  (right).  The  size  and  color  scheme  for  different 
atoms  are  the  same  as  in  Fig.  1.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Table  3 

Lattice  constant  (in  A),  unit  cell  volume  V (in  A3),  magnetic  moment  of  Fe  atom 
MOM  (in  /i b),  and  total  energy  difference  A E  (in  eV  per  unit  cell)  of  Li3FeP04 
and  Li3FeP04 


a 

b 

c 

V 

MOM 

A  E 

L^FePCU 

Non-spin 

12.884 

7.605 

5.364 

525.21 

- 

Reference 

FM 

13.015 

7.593 

5.421 

435.74 

0.98 

-0.54 

AFM 

12.995 

7.657 

5.468 

544.01 

2.66 

-2.25 

LLFePCL 

AFM 

13.372 

7.646 

5.584 

570.90 

2.52 

- 

A E  is  calculated  with  non-spin-polarized  total  energy  as  reference. 


nized  compactly,  while  the  connection  between  each  chain  is 
quite  loose.  To  form  this  structure,  large  volume  expansion  also 
required  (from  Li2FeP04  to  Li3FeP04,  the  unit  cell  volume 
increased  25.4%),  as  shown  in  Tables  2  and  3.  With  such  a 
structure,  lithium  atoms  have  very  good  mobility,  which  will 
be  discussed  in  details  in  the  following  section.  Li  atoms  are 
regularly  distributed  along  the  chained  structure,  stabilized  by 
the  electrostatic  attraction  from  the  PO3  groups.  The  PO3  groups 
are  quite  separated,  which  provide  even  more  space  for  lithium 
storage. 

Recent  experiments  shown  that  lithiation  of  transition  metal 
oxides  leads  to  formation  of  M/L^O  nanocomposite,  acco- 
manied  by  large  volume  expansions,  and  after  lithium  is 
electrochemically  removed  out  from  the  nanocomposite,  the 
remaining  of  the  nanocomposite  is  still  in  nano-sized  metal  oxide 
particles  and  is  used  for  subsequent  cycles  [21-23] .  In  the  case  of 
lithiation  of  LiFePCL  compound,  our  results  show  that  a  quasi- 
one-dimensional  FeO-PC>3  chained  structure  is  formed,  which 
also  provides  good  media  for  lithium  storage.  To  verify  this, 
we  added  more  lithium  in  the  unit  cell  of  this  structure,  and 
the  chemical  composition  becomes  L^FePCL.  As  expected,  the 


Fig.  3.  Schematic  view  of  the  optimized  atomic  structure  of  LLFePCL.  The  size 
and  color  scheme  for  different  atoms  are  the  same  as  in  Fig.  1 .  (For  interpretation 
of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web 
version  of  the  article.) 


chained  structure  is  not  changed  in  Lf^FePCL,  with  16  Li  atoms 
in  the  unit  cell,  as  shown  in  Fig.  3.  These  Li  atoms  are  dis¬ 
tributed  around  the  chained  structure,  which  ensures  that  lattice 
constants  do  not  change  much  upon  Li  intercalation.  As  it  is 
shown  in  Table  3,  the  unit  cell  volume  increases  only  4.9%  from 
Li3FeP04  chained  structure  to  LL^FePCL. 

4.  Electronic  structures  and  charge  distributions  in 
Lii+A:FeP04  (1  <  x  <  3) 

The  electronic  states  of  Lii+xFeP04  around  the  Fermi  level 
are  originated  from  the  Fe  3d  state,  which  degenerate  due  to  two 
types  of  interactions:  the  crystal  field  and  the  exchange  inter¬ 
action.  The  energy  gaps  split  by  these  interactions  are  denoted 
as  Ucf  (for  crystal  field)  and  Ux  (for  exchange  interaction). 
The  split  Ux  separates  the  spin-up  and  spin-down  energy  levels, 
while  the  Uqf  separates  the  five  orbits  within  each  spin  channel. 
In  the  case  of  LiFePCU,  as  reported  in  Ref.  [19],  the  band  gap 
opens  within  the  energy  ranges  of  the  minority  spin  channel, 
due  to  the  local  quasi-octahedral  crystal  field.  As  GGA  calcu¬ 
lation  substantially  underestimates  the  crystal  field  interactions, 
the  calculated  band  gap  is  also  underestimated.  Later,  GGA  +  U 
calculations  improved  the  accuracy  of  the  calculated  band  gap 
[24] .  All  the  results  in  the  present  work  are  obtained  from  GGA 
calculations.  Fortunately,  although  GGA  always  underestimate 
the  band  gap,  the  occupation  states  of  the  energy  levels  are  well 
reproduced. 

The  oxidation  states  of  Fe  atoms  are  +2  in  LiFePCL,  with 
five  electrons  in  the  majority  spin  channel  and  one  electron  in 
the  t2g  part  of  minority  spin  channel.  However,  in  the  case  of 
Li2FePC>4,  the  oxidation  states  of  Fe  atoms  become  +1.  There 
are  seven  electrons  in  the  Fe  3d  orbit,  with  five  filling  in  the 
majority  spin  channel  and  two  in  the  minority  spin  channel.  As 
shown  in  Fig.  lc,  Fe  atoms  are  coordinated  with  three  O  atoms 
within  the  XT-plane.  The  planar  crystal  field  further  degenerates 
the  Fe  3d  orbital  into  five  separate  parts,  z 2  and  x2-y 2  are  the 
lowest  in  energy  in  the  minority  spin  channel  and  become  filled, 
as  shown  in  Fig.  4.  The  results  show  that  the  degeneration  of 


Fig.  4.  Electronic  density  of  states  projected  on  Fe  3d  orbits  in  Li4aLi4CFeP04. 
The  majority  and  minority  spin  channels  are  presented  in  the  upper  and  lower 
panel,  respectively.  The  Fermi  level  is  aligned  to  0. 
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Fig.  5.  Electronic  density  of  states  projected  on  Fe  3d  orbits  in  Li3FeP04.  Plot 
details  are  the  same  as  in  Fig.  4. 

the  Fe  3d  states  under  the  planar  crystal  field  in  the  case  of 
Li^Li^FePCU  is  completely  different  from  that  of  under  cubic 
crystal  field  as  in  the  case  of  LiFePCL.  As  more  electrons  take 
the  minority  spin  channel,  the  calculated  magnetic  moment  of 
Fe  atom  in  Li4flLi4cFeP04  is  3.03,  which  is  smaller  than  that  of 
Fe  atom  in  LiFePCL,  as  listed  in  Tables  1  and  2. 

Presented  in  Fig.  5  is  the  electron  DOS  projected  on  Fe  3d 
orbits  in  Li3FeP04.  From  the  crystal  structure  of  Li3FeP04, 
Fe  atom  is  connected  with  one  P  and  two  oxygen  atoms.  The 
arrangement  of  Fe,  O,  and  P  atoms  is  not  symmetric  to  x-9  y-  or 
z-axis  in  the  unit  cell.  The  DOS  projected  on  Fe  3d  orbits  cannot 
be  clearly  separated.  However,  we  can  see  that  the  occupation  of 
the  minority  spin  channel  is  quite  low,  which  indicates  that  Fe 
atoms  do  not  collect  much  charge  after  Li  atoms  are  intercalated 
into  the  Li2FeP04  lattice  and  become  Li3FeP04. 

To  understand  the  electronic  structure  of  Li2FeP04  and 
Li3FeP04,  it  is  necessary  to  know  the  charge  distributions  in 
the  lattice.  That  is  to  say,  we  need  to  know  when  lithium  atoms 
are  added  into  the  lattice,  how  the  charge  is  redistributed  within 
the  lattice.  For  this  purpose,  we  calculated  the  differential  charge 
density  paiff,  which  is  defined  as  the  difference  between  charge 


density  of  ptLi^FePCL]  and  p[LinFePC>4].  Here  p[LinFePC>4] 
denotes  the  charge  density  after  removing  all  Li  atoms  from 
LixFeP04  without  relaxing  the  atomic  positions.  By  analyzing 
such  a  plot,  it  can  be  seen  clearly  where  the  charge  from  the  Li 
atoms  goes,  as  shown  in  Fig.  6.  In  all  cases,  Li  atoms  are  com¬ 
pletely  ionized  and  in  the  form  of  Li+.  In  LiFePCL,  the  charge 
from  Li  atoms  is  mainly  collected  by  Fe  atoms,  and  thus  the 
oxidation  state  of  Fe  becomes  +2  (see  Fig.  6a).  For  the  case  of 
Lh^Lh^FePCL,  however,  Fe  atoms  do  not  collect  much  share; 
the  majority  part  of  the  charge  from  Li  atoms  is  collected  by  O 
atoms  that  do  not  belong  to  members  of  the  ring  (see  the  crystal 
structure  in  Fig.  lb).  Situation  for  the  case  of  LL^Li^FePCL 
is  also  different,  in  which  the  charge  is  mainly  collected  by  Fe 
atoms,  and  filled  into  the  z 2  and  x2-y2  orbits.  As  shown  in  Fig.  6d 
for  the  case  of  L^FePCL,  both  Fe  and  O  atoms  collected  some 
charge.  Those  O  atoms  that  bonded  with  P  atoms  get  even  more 
share  than  those  connected  with  Fe  atoms.  Charge  redistribu¬ 
tion  is  accompanied  by  the  structural  changes,  during  the  Li 
intercalation  process. 

5.  Migration  of  Li  ions  in  Li3FeP04 

It  is  reported  that  lithium  diffusions  in  LiFePCL  is  one¬ 
dimensional,  which  is  partly  responsible  for  the  poor  rate 
performance  of  the  material  as  cathode  for  lithium  ion  batteries 
[5,7].  However,  the  situation  is  quite  different  for  LiFePCL  as 
anode  material.  We  show  in  this  part  that  migration  of  Li  ions 
in  the  chained  structure  L^FePCL  is  three-dimensional.  For  the 
convenience  of  comparison,  we  first  calculated  the  energy  barrier 
of  lithium  diffusion  in  the  one-dimensional  pathway  in  LiFePCL . 
The  energy  barrier  calculated  with  NEB  method  is  0.29  eV  in 
the  present  work,  which  is  in  good  agreement  with  Ref.  [5]. 

Fig.  7  presents  the  migration  pathway  and  energy  barriers 
for  lithium  diffusion  in  the  chained  L^FePCL.  Five  Li  sites  are 
marked  by  numbers  from  1  to  5  in  the  unit  cell,  as  shown  on 
the  left  side  of  Fig.  7.  There  are  four  diffusion  pathways  are 
simulated:  from  site  1  to  site  2  (denote  as  1-2  and  see  Fig.  7a), 
site  1  to  site  4  (denote  as  1-4  and  see  Fig.  7b),  site  1  to  site  3 
(denote  as  1-3  and  see  Fig.  7c)  and  site  4  to  site  5  (denote  as 


Fig.  6.  Differential  charge  density  pdiff  =/0[LiJCFePO4]  —  pfLinFePCL]  for  (a)  LiFePCL,  (b)  Liz^Li^FePCL,  (c)  Li4aLi4CFeP04,  and  (d)  L^FePC^.  The  color  and 
size  scheme  for  different  atoms  are  the  same  as  in  Fig.  1,  the  isosurfaces  of  the  differential  charge  density  are  in  black  and  transparent.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 
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Fig.  7.  Migration  energy  barrier  for  Li  ions  diffusion  in  chained  Li3FeP04  along  different  migration  pathway.  The  numbers  1-5  refer  to  five  Li  sites  in  the  lattice. 


4-5  and  see  Fig.  Id).  Migration  pathways  of  1-2  and  4-5  are 
migration  of  Li  ions  between  different  chains,  while  migration 
pathways  of  1-3  and  1-4  are  migration  of  Li  ions  within  a  chain. 
Although  the  calculated  energy  barriers  (0.49,  0.46,  0.44,  and 
0.72  eV  for  1-2, 1-4, 1-3,  and  4-5,  respectively)  are  a  little  larger 
than  the  energy  barrier  (0.29  eV)  of  lithium  ions  diffusion  along 
the  one-dimensional  pathway  in  LiFePCL,  they  do  not  differ  very 
much  from  each  other.  This  shows  that  lithium  ions  diffusion  in 
L^FePCL  is  three-dimensional.  Furthermore,  lithium  diffusion 
from  one  chain  to  another  is  more  difficult  than  migration  within 
a  chain,  as  the  energy  barrier  is  a  little  larger  for  lithium  migration 
from  one  chain  to  another. 

6.  Summary  and  conclusions 

In  summary,  Lii+xFeP04  (1  <x<3)  as  anode  material  for 
lithium  ion  batteries  has  been  investigated  using  ab  initio  cal¬ 
culations.  Large  structure  changes  accompanying  large  volume 
expansions  are  observed  when  first  two  Mole  lithium  atoms 
are  intercalated  into  the  LiFePCL  lattice.  Then  the  final  prod¬ 
uct  of  L^FePCL  possesses  a  chained  structure,  which  is  quite 
stable  and  in  favor  of  lithium  storage  and  migration.  Upon  addi¬ 
tional  lithium  atoms  intercalated  into  the  L^FePCL  chained 
structure,  the  chained  structure  keeps  unchanged  and  the  cell 
volume  expansion  is  very  small.  Within  the  chained  structure, 
lithium  diffusions  are  in  a  three-dimensional  pathway,  which  is 
different  to  the  one-dimensional  pathway  in  LiFePCL  as  cathode 
material. 
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